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ABSTRACT 

Context. The current gravitational wave detectors have reached their operational sensitivity and are nearing detection of compact 
object binaries. In the coming years, we expect that the Advanced LIGO/VIRGO will start taking data. At the same time, there are 
plans for third generation ground-based detectors such as the Einstein Telescope, and space detectors such as DECIGO. 
Aims. We discuss the eccentricity distribution of inspiral compact object binaries during they inspiral phase. We analyze the expected 
distributions of eccentricities at three frequencies that are characteristic of three future detectors: Advanced LIGO/VIRGO (30 Hz), 
Einstein Telescope (3 Hz), and DECIGO (0.3 Hz). 

Methods. We use the Star/Track binary population code to investigate the properties of the population of compact binaries in 
formation. We evolve their orbits until the point that they enter a given detector sensitivity window and analyze the eccentricity 
distribution at that time. 

Results. We find that the eccentricities of BH-BH and BH-NS binaries are quite small when entering the Advanced LIGO/VIRGO 
detector window for all considered models of binary evolution. Even in the case of the DECIGO detector, the typical eccentricities 
of BH-BH binaries are below 10~ 4 , and the BH-NS eccentricities are smaller than 10~ 3 . Some fraction of NS-NS binaries may have 
significant eccentricities. Within the range of considered models, we found that a fraction of between 0.2% and 2% NS-NS binaries 
will have an eccentricity above 0.01 for the Advanced LIGO/VIRGO detectors. For the ET detector, this fraction is between 0.4% and 
4%, and for the DECIGO detector it lies between 2% and 27%. 
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1. Introduction 

As the inte rferometric gravitational wave detectors LIGO 
and VIRGO (lHarrv & the LIGO Scientific Collaboration! l2010t 
lAcernese et al.l I2006I) reach their design sensitivities, the first 
detection of gravitational waves has become more imminent. 
Both detectors will undergo serious improvements to increase 
their sensitivity dSrnifh & LIGO Scientific Collaboration! [2009 ; 



i tv _ 
ISpallicci et ai] |2005) 

It is therefore important to investigate the properties of the 
primary candidate sources for detection, namely compact ob- 
ject binaries. There have been a number of papers dealing 
with several properties of the population of compact binaries 
(e.g.. iNelemans & van den Heuvel feOOU IVoss & Taurisl [2003; 
De P onder & Vanbeverenl l2004t ISipior & Sigurdssonl £ 002: 
Pfahl et all l2005b iDewi et alj I2002L l2005t iBogqmazoy et all 
2007t iKiel et alj|2010b . In~particular. lAbadie et al l d2010h pre- 



sente d their est imated detectio n rates . In addition, the mass spec- 
trum dGondek-Rosinska et al. 2007), a nd even spin properties 
dSchnittmanll2004t iMandel & O'Shaughnessvll2010l) have been 
studied. In this paper, we present yet another aspect of the merg- 
ing compact object binary population: the distribution of the ec- 
centricity. 

From radio observations, we currently know of only six 
compact object binaries with merger timescales shorter than a 
Hubble time, all of them NS-NS (neutron star - neutron star) 
systems, and we know of no BH-NS (black hole - neutron star) 
nor BH-BH (black hole - black hole) system. The known NS-NS 
systems are listed along with their orbital parameters in Table [1] 



The observed NS-NS binaries have merger times 
T mer g k, 100 Myr. However, a significant fraction of the 
population of the merging NS-NS may originate in the so-called 
ultra-compact NS-NS binaries, wh ich have much shorter m erger 
timescales T merg < 100 Myr (Bel czvnski et al.l l2002al) . In 
addition, the small number of known pulsars are indicative of 
there being a significant fraction of very eccentric binaries. 

Since no BH-NS nor BH-BH binaries are known, we can 
only rely on evolutionary considerations when estimating their 
number and properties. It has been found that their number de- 
pends very strongly on the outcome of the common envelope 
phase when the secondary is on the Hertzsprung gap. This phase 
will very likely end up as a merger and the formation of a 
Thorne-Zytkow object. However, it has been demonstrated that 
in a low-metallicity environment the common envelope mergers 
may be (to so me extent) avoided and the BH-BH formation is 
very effective dBelczvnski et al.ll20 10b). 

The eccentricity of a compact object binary may potentially 
be derived by analyzing the inspiral signal, provided that the 
eccentricity is significant. In this paper, we investigate the ec- 
centricity distributions in the frequency band of the currently 
working and future detectors of gravitational waves. For the cur- 
rently working detectors (LIGO and VIRGO), we assume that 
the sensitivity of the detectors will allow us to measure the 
signal for the frequencies starting at 30 Hz. This may not be 
accurate for the current state of these instruments, but it does 
accurately represent the predicted sensitivity of the Advanced 
LIGO/VIRGO detectors . We consider two fu ture detectors: 
the Einstein Telescope dVan Den Broeckll2010h and DECIGO 
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(iKawam ura 200f| ISeto et al.ll200lt) . For the Einstein Telescope, 
we assume that binaries shall be detectable from 3 Hz, and for 
DECIGO we assume that the lowest frequency detectable is 
0.3 Hz. In all cases, these have to be treated as indicative num- 
bers that roughly describe these instruments. 

In section 2, we describe the model used to investigate the 
population of compact object binaries. Section 3 presents the re- 
sults for the current and future gravitational wave detectors. In 
section 4, we summarize and discuss the results. 



set of models in which the magnitude of the NS kicks is lower 
by a factor of 2, to cr = 132.5 km/s, as some observations and 
empirically based arguments seem to indica te that natal kicks in 
close binaries are lo wer than for single stars dDessart et al.f 20Q6: 
Kitau ra et al. 2006). The BH kicks are decreased in the similar 
fashion as in models with the full NS kicks. The standard value 
of cr parameter is denoted by K and the smaller value by k. The 
detailed list of models considered in this paper is presented in 
Table|2] Model AZK is a s t andard set of parameters described in 
detail bv lBelczvnski et alJ d2002bl) . 



2. The model 

2.1. Compact object binary population model 

To model the population of compact obj ect binaries, we used th e 
Star Track population synthesis code (Belczvnski et al. 2002b). 
It perform a suite of Monte Carlo simulations of the stellar 
evolution of stars in environments of two typ ical metallicities: 
Z = Zp = 0.02 and Z = 10% Z = 0.002 (e.g.. lBelczvnski et al l 
2010b). In these ca lculations, we employed the recent estimates 
of mass loss rates dBelczynski et alJfeOlOal) . We calculate a pop- 
ulation of 2 million massive binary stars, tracking the ensuing 
formation of relativistic binary compact objects: double neu- 
tron stars (NS-NS), double black hole binaries (BH-BH), and 
mixed systems (BH-NS). Our modeling utilizes updated stel- 
lar and b inary physics, including results from supernova sim- 
ulati ons (|Frver & Kalogeral [2001b and compact object forma- 
tion (iTimmes et al.lll996l) . incorporating elaborate mechanisms 
for treating stellar inter actions such as mass transfer episodes 
dBelczynski et al.l 120081) or tidal synchronization and circular- 
ization dHullll98ll) . We pl ace special em phasis on the common 
envelope evolution phase ( Webbinkl l 19841) . which is crucial for 
close double compact object formation because the attendant 
mass transfer permits an efficient hardening of the binary. This 
orbital contraction can be sufficiently efficient to cause the indi- 
vidual stars in the binary to coalesce and form a single highly 
rotating object, thereby preventing additional binary evolution 
and the formation of a double compact object. Because of signif- 
icant radial expansion, stars crossing the Hertzsprung gap (HG) 
very frequently initiate a common envelope phase. HG stars do 
not have a clear entrop y jump at the core-envelope transition 
dlvanova & Taa m 2004); if such a star overflows its Roche lobe 
and initiates a common e nvelope phase, the inspira l is expected 
to lead to a coalescence dTaam & San dauist 2000). In particu- 
lar, it has been estimated that for a solar metallicity environment 
(e.g., our Galaxy), properly accounting for the HG gap may lead 
to a reduction in the merger rates o f BH-BH binaries by ~ 2 - 3 
orders of magnitude (Belczv nski et al.ll2007l) . In contrast, in a 
low metallicity environment th is suppression is much le ss se- 
vere (~ 1 order of magnitude; iBelczvnski et al.l (1201 Obi) ). The 
details of the common envelope phase are not yet fully under- 
stood, thus in what follows we consider two set of models, one 
that does not take into account the suppression (optimistic mod- 
els: marked with A), and another that assumes the maximum 
suppression (pessimistic models: marked with B). Solar metal- 
licity and 10% of solar metallicity are labeled with Z and z, re- 
spectively. In the case of NSs, we adopt nata l kick distributions 
from observations of single Galactic pulsars dHobbs et al.ll2005h 
with cr = 265 km/s. However, for BHs we draw kicks from the 
same distribution (but at a lower magnitude), which is inverse 
proportio nal to the amount of fall back expected at BH forma- 
tion (e.g.. iFrver & Kalogerall200ll) . In particular, for most mas- 
sive BHs that form with the full fall back (direct BH formation), 
the amount of natal kick is zero. In addition, we test one more 



Table 2: The list of models of stellar evolution used in the paper. 



Model 


Metallicity 


cr [kms '] 


HG 


AZK 


Zg 


265.0 


+ 


BZK 


/ 


265.0 




AZk 


z 


132.5 


+ 


BZk 


Z Q 


132.5 




AzK 


10% Z Q 


265.0 


+ 


BzK 


10% Z Q 


265.0 




Azk 


10% Z G 


132.5 


+ 


Bzk 


10% Z Q 


132.5 





2.2. Evolution of orbits 

The evolution of the orbit of compact object binary under 
the influence of gravitati o nal rad i ation had been calculated by 
iPeters & Mathews! (119631) : iPetersI d!964l) . In the quadrupole ap- 
proximation, the orbit decays as 



da B 
at a 3 



¥(«) = 



1 +73/24e 2 + 37/96e 4 
(1 -e 2 ) 1 ' 2 



(1) 



where a is the great semi-axis, e is the eccentricity of binary, 
M\ is the mass of the first component, M2 is the mass of second 
component, and 
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5 c 5 M x +M 2 
While the eccentricity decays as 



de 19 p 



dt 



12, 



(1 + 121/304e 2 )e 
(1-e 2 ) 5 / 2 ' 



(2) 



(3) 



Using the above formulae we can express the fundamental 
gravitational wave frequency as a function of the eccentricity 



fcw(e) 



2 (1 -e 2 f 12 121 
- [1+ 304* 



,18/19 



-1305/2299 3/2 
C () ' 



(4) 



where c = {ef l9 [\ + 5^e ] 1305/2299 )(l - e 2 Y\ P is the initial 
orbital period, and fcw(e) is the first non-zero harmonic. The 
gravitational wave frequency is twice the orbital frequency, i.e., 

few = 2forb = p~ • 

We present the evolution of eccentricity as a function of 
gravitational wave frequency in Figure Q] for a binary neutron 
star with components of equal masses of 1.4 M Q . The initial fre- 
quency corresponds to a semi-major axis such that the merger 
time is set to be T merg = 10 4 Myr. Figure Q] contains several dif- 
ferent cases of evolution in the plane stretched by eccentricity 
and gravitational wave frequency. 
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Name 


Port m 


Present e 


Tmerge [Gyr] 


e at 0.3 Hz 


e at 3 Hz 


e at 30 Hz 


Ref. 


J0737-3039A/B 


2.454 


0.088 


0.085 


4.5 x 10" 5 


4 X 10-" 


3.5 x 10-' 


Burgav et al. (2003) 


B2127+11C 


8.05 


0.681 


0.2 


2.9 x 10~ 4 


2.6 x 10" 5 


2.3 x 10- 6 


Anderson et al. (1990) 


J 1906+0746 


3.98 


0.085 


0.3 


2.6 x 10- 5 


2.3 x 10~ 6 


2 x 10~ 7 


Lorimer et al. (2006) 


B1913+16 


7.752 


0.617 


0.3 


2.2 x 10~ 4 


1.9 x 10~ 5 


1.7 x 10~ 6 


Weisbers & Tavlor (2005) 


J1756-2251 


7.67 


0.181 


1.7 


2.6 x 10" 5 


2.5 x 10" 6 


2.2 x 10~ 7 


Faulkner et al. (2005) 


B1534+12 (=J1537+1155) 


10.098 


0.274 


2.7 


3.6 x 10" 5 


3.2 x 10~ 6 


2.8 x 10~ 7 


Wolszczan(1991) 




f GW f GW f GW f GW 



Fig. 2: The properties of the population of double neutron stars obtained using the StarTrack code. The plot shows only the binaries 
that will merge within the Hubble time. Solid lines correspond to evolutionary tracks for initial gravitational waves frequencies from 
fo = 10~ 8 Hz (first line from the left-hand side) to fo = 10 2 Hz (first line from the right-hand side). 
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'gw 

Fig. 1: We present ten cases of eccentricity evolution, starting 
with different values of e from e — 0.1 (first line from the 
bottom) to e — 0.99 (first line from the top). Initial semi- 
major axis is chosen such that a binary will merge within time 
Tmerg = 10 Gyr in each case. 

3. Results 

3.1. Properties of the binaries at formation time 

We start with an initial population created using the StarTrack 
code. We present the properties of the population of compact 



object binaries in Figures |2]-|4] in the space spanned by the ini- 
tial eccentricity and initial gravitational wave frequency, which 
is twice the orbital frequency. Each panel in theses figures cor- 
responds to a different model labeled as listed in Table [2] 

The case of the NS-NS systems is shown in Figure [2] The 
boundary of the region populated by the systems on the left-hand 
side corresponds to the requirement that we only consider bina- 
ries that merge within a Hubble time. The bulk of the binaries 
shown in each panel correspond to those that have undergone 
one CE phase in their evolution. The top row corresponds to the 
models AZK, AZk, AzK, and Azk, in which we allow the bina- 
ries to cross through the common envelope with the donor on the 
Hertzsprung gap, denoted by "+" in Table[2] These binaries may 
undergo a second common envelope phase with a helium star 
companion. At the second CE stage, the orbit is tightened even 
more leading to formation of the stripe in the diagram stretching 
from few ~ 10~ 2 Hz at e m 10~ 2 . In these models, the initial 
distribution in the space of gravitational wave frequency versus 
eccentricity is bimodal. The influence of the value of the kick ve- 
locity has a small impact on the shape of distributions presented 
in Figure|2]as can be seen by comparing the data in plots labeled 
as either K-large kicks or k-small kicks. 

For BH-NS systems, presented in Figure [3] and BH-BH bi- 
naries, in Figure |H we present the results of six out of eight 
models, since in models BZK and BZk, almost no binaries are 
formed in our simulations that involve 2x 10 6 initial binaries. For 
BH-NS and BH-BH binaries, the formation of ultra-compact bi- 
naries is not expected. The formation of NS-NS ultra-compact 
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Fig. 3: The properties of the population of neutron star - black hole systems obtained using the StarTrack code. The plot shows 
only the binaries that will merge within the Hubble time. Solid lines correspond to evolutionary tracks for initial gravitational wave 
frequencies from /o = 10~ 8 Hz (first line from the left-hand side) to /o = 10 2 Hz (first line from the right-hand side). 



systems in very close orbits is the consequence of the final CE 
episode, which is initiated by a low-mass helium (2-4 M ) star 
and its NS companion (1.4 M ). Since the donor is about twice 
as massive as its companion, the CE phase is initiated by the 
non-stable mass transfer and the orbit significantly decreased in 
size. For more massive BH-BH/BH-NS binaries, helium stars 
are on average more massive (M > 3 - 4 M ) and do not expand 
(so no CE phase), and even if a low mass helium star forms, 
then its companion is a BH (M > 3 M ), so most likely in- 
stead of CE the RLOF is stable and does not lead to orbital de- 
cay (mass ratio ~ 1). Very few systems (e.g., models AzK or 
Azk) produce ultra-compact BH-NS/BH-BH binaries for very 
special cases of binary evolution. In the case of BH-BH bina- 
ries, shown in Figure [4] we present only six models, since mod- 
els BZ K and BZk do not lead to the formation of BH-BH bi- 
naries dBelczvnski et alj 120071) . In all models, there is an en- 
hanced density of systems formed with e « 0.1 at approximately 
lCT 5 Hz < few < lCT 4 Hz. In these systems, the second black 
hole has formed via direct collapse. When treating the direct col- 
lapse, we assume that 10% of the mass escapes in the form of 
neutrinos and possibly gravitational waves. Hence, the gravita- 
tional mass of the BH is 10% lower than the baryon mass of the 
collapsing star. This introduces a small eccentricity a; 0.1 since 
the systems were circularized in the mass transfer prior to the 
collapse and the formation of the second BH. 

3.2. Eccentricity when binary enters detector band 

For the detection of gravitational waves, it is important to know 
the eccentricity of a binary at the time it enters the sensitivity 
window of the detector. We consider three cases that correspond 
approximately to three types of detectors. In Figure [5] we show 
the eccentricity distributions at 0.3 Hz (bottom horizontal axis) 



and 3 Hz (top horizontal axis) corresponding approximately to 
the ET and DECIGO detectors. The results for the Advanced 
LIGO/VIRGO can be easily obtained by rescaling the horizontal 
axis. 

The shape of the eccentricity distributions at the moment that 
the binary enters the given detector band follows from the cor- 
responding initial distribution. However, one must note that for 
each type of binary there is a different natural timescale and fre- 
quency, because of the different mass scales of each binary. 

We present the results for the DECIGO detector and add ap- 
propriate numbers for the ET in parentheses. For the NS-NS bi- 
naries shown in the top panel of Figure [5] the distribution is ei- 
ther centered ones 10~ 4 (ET: 10~ 5 ) for the models BZK, BZk, 
BzK, and Bzk, where we do not allow the formation of ultra- 
compact binaries in a second CE phase. The remaining mod- 
els AZK, AZk, AzK, and Azk contain another component cen- 
tered roughly at e m 10~ 4 (ET: 10~ 3 ). This additional compo- 
nent represent the ultra-compact binaries that have experienced 
two episodes of mass transfer in their evolutionary history and 
were already very tight at the second supernova explosion. The 
mixed BH-NS binaries, shown in the middle panel of Figure [5] 
exhibit a distribution of eccentricity centered ate * 3x 10~ 5 (ET: 
3 x 10~ 6 ), while the eccentricity BH-BH binaries, shown in the 
bottom panel of Figure [5] lie between e as xlO~ 6 (ET: 10~ 7 ) and 
e a xl0~ 4 (ET: 10~ 5 ) . 

For the Advanced LIGO/VIRGO detectors where we assume 
that the low frequency boundary lies at m 30 Hz, the eccentric- 
ities are even smaller. It follows from equation [4] that the dis- 
tributions are shifted by a factor of 10 -19 ^ 18 for each factor of 
ten in frequency. Thus, the values of eccentricity in the case 
of Advanced LIGO/VIRGO type detectors are consistent with 
e — and we can safely assume that all BH-NS and BH-BH 
binaries are circular without any loss of sensitivity. 
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Fig. 4: The properties of the population of double black holes obtained using the StarTrack code. The plot shows only the binaries 
that will merge within the Hubble time. Solid lines correspond to evolutionary tracks for initial gravitational waves frequencies from 
fo = 10~ 8 Hz (first line from the left-hand side) to fo = 10 2 Hz (first line from the right-hand side). 



In table [3] we present the fraction of binaries with eccentric- 
ities above 0.01 at the time of entering the detector band, to help 
quantify the extent of the large eccentricity tails of the distribu- 
tions presented in Figure [5] T his fraction does not reflect the de- 
tectability of eccentricity (Shapir o Key & Cornis h 2010),which 



for realistic distributions of binaries will be discussed in a forth- 
coming paper. 

4. Summary 

We have presented the eccentricity distributions of compact ob- 
ject binaries at three frequencies immediately before merger. 
The properties of the compact object binaries have been calcu- 
lated using the StarTrack population synthesis code. We have 
found that the eccentricity distributions of the compact object 
binaries do not depend strongly on the assumed model of binary 
evolution. Any dependence has been found to be the strongest 
for binary neutron stars, whose distributions may be either single 
or double peaked. The extra peak corresponds to ultra-compact 
NS-NS binaries that have undergone an additional CE phase im- 
mediately before forming the second NS. 

To make the results easier to use in the simulations, we have 
fitted the resulting distributions of eccentricity with a single log- 
normal distribution in the case of BH-BH and BH-NS binaries 

\2\ 

(5) 



/(*) = 



1 



V2^ 



exp - 



2cr 2 



where x = log <?, fi is the mean, and cr is the variance. 

For the NS-NS eccentricities, we used a sum of two log- 
normal distributions with two weights, since the distribution is 
double peaked, given by 

(x-^f) , (1-w) { (x-^ 



/(*) 



V2^ 



exp 



2<x 2 



(1-w) 
o"2 y/2n 



exp 



2a\ 



(6) 



Table 3: Fraction of the compact binaries with eccentricity 
greater than 10 -2 . Top table corresponds to double neutron stars 
(NS-NS), middle to the mixed systems (BH-NS), and the bot- 
tom to the binary black holes (BH-BH). We present the fraction 
of binaries at the moment of entering different frequency bands 
(30 Hz, 3 Hz, and 0.3 Hz). In brackets, we include the number 
of these systems in the simulation N. We only listed results for 
models that are non-zero. The number of digits shown is for for- 
matting only, and the relative sampling error is N~ 1 ^ 2 . 

NS-NS = 





30 Hz 


3 Hz 


0.3 Hz 


AZK 


0.60% 


(51) 


1.32% 


(112) 


11.13% 


(945) 


BZK 


1.27% 


(36) 


2.33% 


(66) 


6.52% 


(185) 


AZk 


0.16% 


(27) 


0.38% 


(64) 


10.37% 


(1732) 


BZk 


0.30% 


(15) 


0.75% 


(37) 


2.22% 


(110) 


AzK 


0.29% 


(25) 


0.96% 


(83) 


21.74% 


(1880) 


BzK 


1.87% 


(13) 


4.02% 


(28) 


9.33% 


(65) 


Azk 


0.26% 


(37) 


0.57% 


(81) 


26.91% 


(3799) 


Bzk 


1.74% 


(21) 


3.31% 


(40) 


7.79% 


(94) 


BH-NS 




30 Hz 


3 Hz 


0.3 Hz 


AZK 


0.29% 


(2) 


0.73% 


(5) 


3.05% 


(21) 


AZk 


0.15% 


(2) 


0.54% 


(7) 


0.61% 


(8) 


AzK 


0.56% 


(14) 


0.96% 


(24) 


3.96% 


(99) 


BzK 


0.68% 


(10) 


1.23% 


(18) 


3.63% 


(53) 


Azk 


0.35% 


(15) 


0.78% 


(34) 


2.81% 


(122) 


Bzk 


0.33% 


(8) 


0.91% 


(22) 


1.53% 


(37) 


BH-BH 




30 Hz 


3 Hz 


0.3 Hz 


AZK 0.3F 


t (1) 


0.62% 


(2) 


1.87% 


(6) 


AzK 0.02% (3) 


0.02% 


(4) 


0.13% 


(23) 


BzK 0.15% (2) 


0.15% 


(2) 


0.46% 


(6) 


Azk 0.00% (1) 


0.01% 


(2) 


0.03% 


(6) 
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Table 4: Parameters of log normal distribution fitted to results of 
model AZK with asymptotic standard errors. 



-5 -4 -3 -2 
log(e) DECIGO 

log(e) ET 
-5 -4 



log(e) ET 



-3 




-5 -4 -3 -2 -1 
log(e) DECIGO 



-2 



I I 

~ ,'/' \ 

</ 1 

■* ' fx 7 \ 

t> /\/u 


1 

K 

! A 


1 1 1 
AZK - 

AZk - 

AzK . 

BzK 

Azk 

Bzk _ 


1,1, 


^\ 

, 1 


* i i- ■ - 



-7 -6 -5 -4 -3 -2 -1 
log(e) DECIGO 

Fig. 5: Distribution of eccentricity for (NS-NS - top panel, BH- 
NS middle panel and BH-BH bottom panel) seen at 0.3 Hz 
(DECIGO -like detectors) and at 3 Hz (ET - like detectors). Solid 
thick line corresponds to standard model (AZK). Dashed and 
dotted lines indicate other models. 



NS-NS 




0"i : 


0.47 ± 0.04 






o- 2 : 


0.39 ± 0.03 






w. 


0.46 ± 0.03 






0.3 Hz 


3 Hz 


30 Hz 


Hi 


-3.83 ± 0.04 


-4.89 ± 0.04 


-5.94 ± 0.04 




-2.32 ± 0.03 


-3.38 ± 0.03 


-4.43 ± 0.03 


BH-NS 




cr: 


0.55 ± 0.03 






0.3 Hz 


3 Hz 


30 Hz 




-4.18 ± 0.04 


-5.23 ± 0.03 


-6.27 ± 0.03 


BH-BH 




cr: 


0.70 ± 0.03 






0.3 Hz 


3 Hz 


30 Hz 


M 


-4.91 ± 0.08 


-5.95 ± 0.03 


-7.06 + 0.03 



the neutron star. Second, the initial kicks at formation of BHs 
are lower than in the case of NS, so the initial eccentricities of 
BH-BH systems are typically lower than in the case of NS-NS 
ones. 

The eccentricities of the mixed BH-NS systems are larger 
than in the case of BH-BH ones. However, the number of sys- 
tems is small enough to ensure that by neglecting the eccentricity 
we do not decrease the sensitivity of Advanced LIGO/VIRGO 
detectors. For the ET-like detector, some eccentric systems may 
be detected. In the case of the DECIGO-like detector, the num- 
ber of systems with eccentricities above 0.01 lies between 3% 
and 4%. 

The eccentricity of NS-NS systems are larger than those of 
binaries containing BHs. Given a much larger expected detection 
rate for ET, this means that there should be a significant number 
of NS-NS binaries with detectable eccentricities. Finally, in the 
case of DECIGO a fraction of between 2% and 27% of the NS- 
NS binaries have eccentricities above 0.01. Moreover, the shape 
of the eccentricity distribution of NS-NS binaries will depend 
on the existence of an evolutionary scenario leading to the for- 
mation of ultra-compact binaries. Thus, the measurement of the 
eccentricity distribution is an interesting tool for probing the de- 
tails of NS-NS formation scenarios. 
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where x = log e, /j\ is the mean of the first peak, /j.2 is the mean 
of the second peak, cr\ is the variance of the first distribution, a~2 
is the variance of the second distribution, and w is the weight. 

We used the Marquardt-Levenberg algorithm to find the pa- 
rameters, and estimate the asymptotic standard error of each of 
them. The results of the fits are shown in Table |4] The widths 
of the distributions cr are the same for each frequency band and 
only the centroids move. 

The eccentricity of the BH-BH binaries in all three cases is 
negligible. This is due to two factors. First, at the frequencies of 
interest the BH-BH systems are much closer to coalescence than 
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